Hall of Mirrors Scattering from an Impurity in a Quantum Wire 



J. Y. Vaishnav/'B A. Itsara, 2 and E. J. Heller 1 ' 3 

1 Harvard University Department of Physics, Cambridge MA 02138 
2 Harvard College, Cambridge MA 02138 
'Harvard University Department of Chemistry and Chemical Biology, Cambridge MA 02138 

This paper develops a scattering theory to examine how point impurities affect transport through 
quantum wires. While some of our new results apply specifically to hard-walled wires, others-for 
example, an effective optical theorem for two-dimensional waveguides-are more general. We apply 
the method of images to the hard-walled guide, explicitly showing how scattering from an impurity 
affects the wire's conductance. We express the effective cross section of a confined scatterer entirely 
in terms of the empty waveguide's Green's function, suggesting a way in which to use semiclassical 
methods to understand transport properties of smooth wires. In addition to predicting some new 
phenomena, our approach provides a simple physical picture for previously observed effects such as 
conductance dips and confinement-induced resonances. 
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I. INTRODUCTION 

Elastic scattering from a point defect in a hard- 
walled, multimode quantum waveguide is a problem 
which has been the subject of experimental^ and 
theoretical 3 ^**!^^*^^ inquiry in both condensed matter 
and atomic physics. In this paper, we revisit the problem 
from an unconventional point of view, reducing it to the 
scattering of a single effective wavefunction off an array 
of images. Our approach enables us to easily understand 
complex transport phenomena which have already been 
observed (such as conductance dips), as well as predict 
other new ones. 

The problem of scattering from pointlike impurities in 
waveguides is a general one, and thus has physical appli- 
cations in several fields. Impurity scattering in quantum 
wires, and the resulting disorder effects on electron trans- 
port, have long been of interest in mesoscopic physics. 
Specifically, the problem arises in the context of two- 
dimensional electron gases (2DEGs) confined by quan- 
tum wire potentials; the impurity can represent a defect 
in the wire. Such quantum waveguides can be fabricated 
at low-temperature Gai-^Al^As interfaces, and are of 
interest due to their potential role in high-frequency and 
quantum devices. As carbon nanotubes can behave as 
few-mode quantum waveguides^ our results are relevant 
to transport in nanotubes with adsorbed impurities-of 
interest in the study of biosensors. The problem of scat- 
tering in confined geometries arises in atom waveguides 
as well, with potential implications for quantum comput- 
ing and atom intcrfcromctry. In both cases, atoms must 
maintain coherence as they pass through the waveguide, 
and the goal is thus to minimize phenomena such as colli- 
sional phase shifts^ Finally, the wave phenomena arising 
in quantum waveguides are in direct mathematical cor- 
respondence with a number of wave phenomena in other 
systems: Water waves resonantly trapped between an 
array of cylinders^ large antenna arrays^ and sound 
waves directed by "Bessel" line arrays in acoustics^ 

Motivated by the above applications, we examine low- 



energy scattering of noninteracting particles from impu- 
rities in an infinite two-dimensional quantum wire, with 
hard walls at y = and y — d. The wire (Fig. ^) con- 
tains a point impurity, which we model as an s wave 
scatterer of effective radius a. Our method differs from 
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Figure 1: Schematic of the quantum wire, with a point impu- 
rity of scattering length a at fo = (xo, yo)- We assume infinite 
leads, and hard walls at y = and y = d. 

previous treatments, such as the renormalized t matrix 
method) 12 i 16 i 17 in that we use an unconventional ap- 
proach which combines the following three ideas: 

1. The method of images. 

2. The realization that a confined scatterer, like a free 
space scatterer, is a rank-one target. We can, there- 
fore, combine the N degenerate transverse modes 
into a basis in which only a single wavefunction 
scatters. 

3. That the effect of scattering from the waveguide 
walls is to renormalize this single scattering wave- 
function, so that it takes on a new effective value 
at the scatterer location. 

Using the method of images allows us to explicitly un- 
derstand how each reflection from the waveguide walls 
affects the conductance. The fact that the confined scat- 
terer is a rank one target enables us to express the trans- 
port properties of the wire entirely in terms of the single 
scattering wavefunction. Finally, the idea of renormaliza- 
tion allows us to fully explain the transport properties of 
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the wire in terms of reflections of this single wavefunction 
from the waveguide walls. 

In Section [H] we review free space multiple scattering, 
anticipating phenomena which will appear in the wire. 
Section IIIII lays out a similar formalism for scattering 
in the waveguide. In Section II VI we present physical 
phenomena which we observed in our model, which we 
interpret in Section^in terms of interference effects. We 
relegate mathematical details to Appendices IAIBI 



A. Single Scatterer 

Consider a single scatterer at r — tq. In our t ma- 
trix representation, an incident wave <j>(r) scatters into a 
wavefunction ip(f) according to the Lippmann-Schwinger 
equation 

= 4>if) + 8(k)(f>(f6)Go(f, r - k) (5) 

where 



II. PRELIMINARIES: REVIEW OF FREE 
SPACE SCATTERING 



In Section IlIII we show that the Lippmann-Schwinger 
equation for scattering from a confined impurity resem- 
bles the free space Lippmann-Schwinger equation, with 
the exception that in the wire, a renormalized incoming 
wavefunction replaces the true one. Surprisingly, even 
for a single confined scatterer, this renormalized wave- 
function includes free space multiple scattering effects. 
We thus begin by briefly reviewing free space scattering, 
emphasizing the three phenomena which will reappear in 
the Lippmann-Schwinger equation for a single scatterer 
in the waveguide. 

Consider N pointlike impurities in two-dimensional 
free space, located at {fi, . . . ,fjv} . Let V(\r — f;|) be the 
potential at f due to the i th scatterer. We seek scatter- 
ing solutions tp(r) to the time-independent Schrodinger 
equation 



N 



2m 



V 2 + Y,V(\r-r l \)~ E 



W) = 0, (1) 



under the assumption of low-energy, s-wave scattering. 
Except as noted, we henceforth use atomic units, where 
h = m = c = 1. Applying the t matrix formalism^ we 
characterize a single scatterer at by its t matrix, 



t = s(k)\r i ){? i \ 



(2) 



where s(k) is a function of the wavenumber, k = 
v2rnE/H. We choose the functional form of s(k) to sim- 
ulate the scatterer of interest, under the constraint that 
s(k) must satisfy the optical theorem 



G (f,fo; k) = 



2m 



- l -H^\k\r-r \ 



(6) 



is the 2D free space advanced Green's function satisfying 



(V 2 + k 2 ) G {r,f ;k) = ^ 2 \t~ r ). (7) 

We shall henceforth assume we are using the advanced 
Green's function and omit the superscript on the Hankcl 
function; we shall also omit the implicit k dependence in 
the Green's function. 
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Single Scattering Wavefunction 

Any incoming wave <fi(r) must be a solution to the free 
space Schrodinger equation. We can therefore express 
4>{r) in the basis of "cylinder harmonics," 



E 

m=0 



J m (k\r- fb|) x 



Jm(4>-4> ) 



(8) 



where m = 0, 1, 2, . . . correspond to s,p,d, . . . waves. An 
s wave scatterer in free space is a rank one perturbation: 
Of all the terms in (JSJ) , only the m = term is nonzero at 
the scatterer. From (jSJ), then, only the s wave scatters, 
acquiring a phase shift; remaining higher partial waves 
pass through the scatterer unperturbed. While trivial in 
free space, these facts will reappear more subtly in the 
wire. 



-2Ims(fc) = |s(fc)| 2 . 
In this paper, we choose 

h 2 Jo(fca) 



s(k) 



-2i- 



which represents a hard disk of scattering length a. 



(3) 



(4) 



B. Multiple Scatterers 

The Lippmann-Schwinger equation describing scatter- 
ing from a collection of N identical point scatterers at 
{fi, . . . ,f N } is 



N 



(9) 
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in which we express the il>i(fi) recursively as 



= fan 



N 

3 = 1 

3¥* 



(10) 



where we have followed Foldy's method^ Comparing 
II Op t with the single-scattcrcr version (JjjJ, the tpi(fi) are 
the effective incoming wavefunctions at each scatterer: 
tpi(fi) is the amplitude incident on the i th scatterer after 
scattering from each of the other scatterers. A crucial 
point is that ipi(fi) excludes the singular self- interaction 
of the scatterer. 

Defining fa = fafi), ipi = ipffi), inverting (|9I10|I yields 



ip = (1 - sG) 



where 



Gi 



G (ri,fj) i ^ j 
i=j 



(11) 



(12) 



excludes the singular term. 
Ulljl is as the Born series 



An alternate expression of 



$ = [1 + sG + (sG) 2 + (sG) 3 + . . .] 



(13) 



The effective wavefunctions thus have a simple interpre- 
tation in terms of interfering paths: The terms in square 
brackets describe amplitude incident at r*i after interac- 
tions with zero, one, two, or three other scatterers re- 
spectively. The series continues infinitely. 

Anticipating phenomena that will reappear in the wire, 
we highlight the following points of this section: 

1 . The effect of multiple scattering is to create a new 
effective incoming wavef unction at each scatterer. 

2. The effective wavefunction at a particular scatterer 
is a sum of waves scattered from all the other scat- 
terers, and excludes the (singular) self-interaction 
of the scatterer. 

3. An s wave scatterer in free space is a rank one tar- 
get. 



III. SCATTERER IN A WIRE 

In this section, we examine how confinement in a wire 
affects the transport properties of a scatterer. Applying 
the method of images, we derive a form of the empty wire 
Green's function. Using this form of the Green's function, 
we show that the sole effect of confinement is to renor- 
malize the incoming wavefunction. The image represen- 
tation allows us to very simply calculate the renormaliza- 
tion coefficient. We define an effective cross section for 
a confined scatterer, and relate it to the conductance of 
the wire-impurity system. We derive an effective optical 



theorem in the wire. Finally, applying this formalism, we 
investigate the behavior of the cross section and conduc- 
tance as functions of various parameters. 



A. Green's Function via the Method of Images 

In order to write down the Lippmann-Schwinger equa- 
tion for an impurity in the wire, we require the Green's 
function G w (f,fo) of the empty wire. The usual spectral 
form of the empty wire Green's function is 



G w {r,r ) 



1 



^ L.("0 

m—l r»x 



Xm(y)xm(yo)e 1 



(14) 



where the Xm{y) are the transverse modes for the partic- 
ular waveguide (see e.g. DattaS). In our hard wire, the 
modes have the form 



Xm(y) = 







fsin(^) 0<y<d 

y<0,y>d. 



(15) 



We note that the spectral Green's function l|14|) is a sum 
over evanescent modes as well as propagating ones. 

Although we could, in principle, proceed using the 
spectral form Ijl4|l of the Green's function, the spectral 
Green's function provides little physical insight into the 
precise mechanism of the scattering. We therefore take 
an alternative approach, and use the method of images to 
derive an alternate, equivalent form of i|14|) for the hard 
wire. 

Although the problem of an impurity in a hard 
wire has been the subject of numerous theoretical 
studies^*^*S*2i2ii£ the authors are aware of only one 
work 3 , which applies the method of images: Ref. 3 treats 
the related problem of a point scatterer in a 3D waveg- 
uide of rectangular cross section, which the authors re- 
duce to the problem of scattering from an impurity in a 
finite, 2D box. The method of images has the significant 
advantage of making multiple scattering effects explicit. 
We shall thus use the approach of Ref. 3, rather than the 
more common spectral methods. Because our problem of 
an infinite waveguide is somewhat simpler than the finite 
box treated in Ref. 9, the role of multiple scattering is 
more transparent. We are therefore able to explore in de- 
tail the considerable effect of multiple scattering on the 
physical properties of the wire, which, although briefly 
mentioned, is not examined in Ref. 3. 

The empty wire Green's function G w (f, r*o) satisfies 
Green's equation inside the wire, and is zero on the wire 
walls: 



{V 2 + k 2 )G w (r,f ) = 25(r~r ) (16) 
G w (x,0) = G w (x,d)=0 (17) 

The confined Green's function differs from the free space 
version due to the necessity of including reflections off the 
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walls when describing the response to a point excitation. 
Note the factor of two in (|lfc>|) , obtained by casting (JJJ in 
atomic units. 

As in electrostatics, we can linearly combine free space 
Green's functions © with point sources in different 
places, so that the linear combination satisfies i|17|) . As 
long as only one of the point sources is inside the wire 
(0 < y < <£), our sum will satisfy inside the wire. 
Outside the wire, we discard the solution. If the point 
source is located at r = (xq, yo), we use the configuration 
of image point sources in Fig. We reflect the scatterer, 
creating a series of images such that the position of the 
n th image is r n = xq£ + [(— l) n yo + 2nd]y. 
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Figure 2: Fhe method of images allows us to reduce the prob- 
lem of a single scatterer in a wire to that of a periodic array 
of scatterers in free space. Fhe solid lines and black scat- 
terer are the actual wire and point source. Fhe gray point 
sources, and the dashed lines, are images. Fhe pluses and mi- 
nuses refer to the sign of the free-space Green's function term 
that the particular point source contributes. Fhe alternating 
signs of these contributions are due to the Dirichlet boundary 
conditions; the wavefunction must cancel on the wire walls. 

This image configuration yields an empty wire Green's 
function of the form 

oo 

G w (r,r )= ]T (~l) n G (f/ n ) (18) 

n— — oo 

The Green's function (|18fl satisfies (|16fl inside the wire, 
and l|17|) on the wire boundaries. Outside the wire we 
implicitly set the Green's function to zero. The Green's 
function satisfies the Dirichlet boundary conditions be- 
cause, as we take infinitely many images, the alternating 
sum of Hankcl functions converges to zero on the walls. 
This series converges extremely slowly, precluding its use 
for numerical purposes (see the discussion in Morse and 
Feshbacbja). 



B. Renormalizing the Incoming Wavefunction 



The Fippmann-Schwinger equation requires knowledge 
of both the Green's function and the incoming wavefunc- 
tion. In Section Hi Al we calculated the Green's function 
for the wire. In this section, we discuss the incident wave- 
function. We show that the effects of confinement resem- 
ble free space multiple scattering: Confinement renor- 
malizes the effective wavefunction at the scatterer from 
its true value. In the image formalism, the multiple scat- 
tering arises due to scattering off images. 

Applying the free space multiple scattering formalism 
(|9I1U|) to the array of images in Fig. |3 we find that the 
Lippmann-Schwinger equation for the wire is 

oo 

V(r) = 0(r) + S J2 G Q (r,rMn) (19) 

i— — oo 

where 

CO 

ipi (n) = 4>(ft) + s G ° ( f i > ?i ft ) ■ ( 2 °) 

J= — oo 

Recall from our discussion of free space multiple scat- 
tering the physical meaning of this recursive Lippmann- 
Schwinger equation: Fhe ipi(ri) defined in (|2U|) are the 

effective incoming wavefunctions at the scatterer in 
the array. In a general multiple scattering problem, we 
would have to stop here and solve numerically. Fhe im- 
age array is, however, a special case: Not only is it peri- 
odic, but the boundary conditions at the wall impose the 
antisymmetry 

<j>(n) = (-i)Wo) 
Mn) = (-i)Vo(fb)- (21) 

Any wavefunction, incident or scatterered, must be zero 
on the walls. Thus, inside the wire, in the y direction, the 
wavefunction is a superposition of modes (|15|) . Replacing 
the wire with an image array requires that we extend 
these modes outside the wire, into the region y < 0, y > 
d. Keeping the Xm(y) continuous at the walls, we can 
extend them so that the wavefunction is either symmetric 
or antisymmetric across the wire walls; either set forms 
a basis. However, we know that in the absence of a wall 
Xm(y) must be continuous also, as we have eliminated the 
hard walls. The symmetric extension does not satisfy this 
condition. We therefore discard the symmetric extension, 
leaving l|21|) as the proper boundary condition. 

Combining (|2U|) and (|21|) , the equation for the wavelets 
thus becomes 

Mn) = (-I) 1 mr a ) + sMro)G r ] (22) 
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where we have defined the renormalization sum 



G r = J2 (-l) 4 Go(rl,r* ) 



(23) 



i= — oo 



noting that 



G r = (-iy i-rfGoirufj)- 



(24) 



l — — oo 



Taking i = in gZ) , we find 



■00 (rb) 



1 - sG r 



(25) 



Substituting I|21I25|I into l|19|) yields the Lippmann- 
Schwinger equation for scattering from an impurity in 
a wire: 



if)(r) = 0(r) + s 



1 - sG r 



1 - sG r 



J2 (-l) i G (r,rl)(26) 

i— — OO 

G w (r,r ). (27) 



Comparing l|27|l to the free space Lippmann-Schwinger 
equation JSJ for a single scatterer, we see that the wire 
simply renormalizes the incoming wavefunction at the 
scatterer to have a new effective value 



(28) 



1 - sG r 

and so our final Lippmann-Schwinger equation is 

0(r) = 0(rO + s4>(f )G w (r,f ). (29) 
For convenience, we define a renormalization factor 



K 



4>(ro) 



(30) 



which is the ratio between the original and renormalized 
incoming wavefunctions at the scatterer. 



In order to cast the renormalized incoming wavefunc- 
tion in terms of interfering paths, we expand perturba- 
tively, as we did in i(T3|) for free space. For shorthand, 
define 

G^^i-iy^Gin^). (31) 

We find 

0(ffe) = (l + sG r + s 2 G 2 r + ...)<b(ro) 



E 

i— — oo 



sG,_ 



E sG '^o x sGj^o 



Z— — 00 J — — 'DO 



0(»*o) 



E 

i— — oo 



(32) 



E 



E ■ 

— — oo J — — oo 



sGi_>j x sGj^o 



where we have used (|24|) . The n th term of the sum 
describes propagation to r*o after scattering from any n— 1 
images, beginning with the first term, which describes 
free propagation to r*o- Alternately, the n th term corre- 
sponds to n scattering events combined with any number 
of reflections off the wall, each of which introduces an ad- 
ditional phase of -1. Fig. ^illustrates the correspondence 
between scattering from an image and reflecting from the 
wall. 




Figure 3: Amplitude scattering from the n th image corre- 
sponds to amplitude reaching a point after n bounces from 
the wall, each of which reverses the sign of the wavefunction 
(phrased semiclassically, reflections introduce a Maslov index 
of -1). Resonances occur when the interference of all paths is 
maximally constructive at the scatterer itself. 



Renormalization of the effective wavefunction is sim- 
ply the manifestation of the identical phenomenon in free 
space (see item^in the list of Section lll B|) . From a semi- 
classical point of view, the incoming wave, given infinite 
time to spread, reflects from each possible combination 
of scatterers before returning to the source. Note that 
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G r is simply the wire's Green's function evaluated at the 
scatterer, with the singular self-interaction of the source 
removed, 



G r = lim [G w (f,f ) 
r— »r 



Go(r,fo)]. 



(33) 



This idea, as expressed in (|3^1) . is the essence of "renor- 
malized t matrix theory. 'U&iSiil The theory, as applied 
to our case, is equivalent to the idea that in free space 
multiple scattering, the effective wavefunction at a scat- 
terer excludes the singular self- interaction (see item |3 of 
Section IIIBfl . We note that although we have derived 
(|33|l via images only for the special case of a hard walled 
guide, the result (usually derived perturbatively) is in 
fact true for an arbitrary guide or confining potential. 16 



Numerics 

While the method of images led quickly to the central 
results (|27I33I) , these results expressed in terms of image 
sums converge so slowly as to be impractical for numerics. 
Although the spectral form l|14|l converges more rapidly 
than the image series Q18[l. its convergence is still not 
uniform, as the evanescent modes include a logarithmic 
singularity. 

In AppcndixlBl we use Kummer's method to accelerate 
the convergence of the Green's function, casting it in the 
form l|B14|) suitable for numerical work. A side benefit 
of applying Kummer's method is that we obtain a more 
rapidly converging expression for G r : In Appendix [51 we 
show that the scries <|23[) for G r (which is a Schlomilch 
series) resums to 



C. Method of Images vs. Spectral Formulation 

We have presented both the image series (|18|) and the 
spectral series (11411 for the Green's function. Because 
the Green's function of a system is unique, these two 
forms must be equivalent. However, the physical rela- 
tionship between the method of images and the spectral 
form is far from obvious. The two forms of the Green's 
function highlight different physical phenomena. For ex- 
ample, reflections from the walls, which appear imme- 
diately in (|18|l . are far less obvious in the spectral form 
(I14JI . Contrastingly, the role of evanescent channels in the 
scattering, while straightforward in the spectral Green's 
function l|14ll . is less transparent in the image expansion 
()18J) of the same Green's function (although even the im- 
age expansion suggests that evanescent channels will be 
present in some form, because Hankel functions are sin- 
gular). In order to understand the connection between 
the image and spectral forms of the Green's function, we 
show their equivalence mathematically in Appendix lAl 

In Appendix A, we use an integral form of the Han- 
kel function to show the equivalence of Ijl8|l and (|I4|I . 
The physical connection between the image and spectral 
formalisms is diffraction. The images in Fig. |3 form a 
periodic array, which is effectively a diffraction grating. 
We note one small difference from a typical diffraction 
situation: Due to the wire boundary conditions, the inci- 
dent modes are superpositions of two plane waves each, 
rather than a single incident plane wave. This difference 
is trivial, due to the linearity of Schrodinger's equation. 
When such modes are incident on the effective lattice of 
images, they diffract: Plane waves (or superpositions of 
plane waves) striking the grating scatter, at large dis- 
tances, into sums of plane waves at the Bragg angles. 
The diffracted spectral orders are evident as the quan- 
tized wavenumbers appearing in the spectral form of the 
Green's function. As the Green's function for each im- 
age scatterer is singular, evanescent waves appear in the 
Fresnel regime, near the lattice. We examine these state- 
ments rigorously in Appendix lAl 



G r — 
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7T 



- Sm {— ) 
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(34) 



where 7 is the Euler-Mascheroni constant. We note that 
this expression is considerably more complex than the 
expression ()23J) derived via the method of images. In 
the image formalism, the expression for the renormal- 
ization constant G r has the intuitive form lll'-'il) . because 
removing the contribution of the source simply involves 
excluding the scatterer itself, while retaining the images. 
In the spectral formalism, because the singularity is not 
explicit, its removal is considerably less transparent. 



D. Scattering Phenomena in Quantum Wires 

In this subsection we define an effective cross section 
and optical theorem for the confined impurity. As we 
have shifted from the image formalism (which applies 
only to the hard wire) to the spectral one, where (|14|) 
applies generally, the results presented in Section UlI Dl 
apply to arbitrary waveguides, not only hard wires. 



1. The S Matrix 

Since the wire is infinite, for calculating transmis- 
sion/reflection at infinity, we need consider only the open 
channels. Without loss of generality (due to translational 
invariance of the wire in the x direction), suppose the 
scatterer is at xq = 0. With incident mode n, normalized 
to unit flux, the wavefunction at large distances from the 
scatterer is 



-Xn{y) 



rn — 1 



ilZs 



An) Am) 
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x Xn(yo)xm(yo)' 



ife£ m) |x| 



-Xm(y) 



(35) 



from which we can read off a scattering matrix 

S= ( ™ ^, ) = ( . "\. " „" i (36) 



where the reflection and transmission coefficients are 

iRs 



1 L™ 



An) Am) 



--Xn(yo)x-m(yo) (37) 



(38) 



One can show algebraically (or verify numerically) that 
rank(R) = rank(T - I) = 1. 



2. Defining a Cross Section in the Wire 

In this section, we define an effective scattering cross 
section for the wire. The usual three-dimensional free 
space differential scattering cross section is 



da 



dN(Q) 
Ni„dQ 



(39) 



where dN(il) is the number of particles per area scattered 
into the solid angle dfl, and N in is the number of incident 
particles per unit area. Expressed in terms of probability 
currents, the above relation becomes 



da J j r r 2 dfl 

d£l jin 



(40) 



where j = Im(^ r *V^ r ) is the probability current, ji n is 
the incoming probability current, and j r is the radial 
probability current. Integrating over solid angles, we find 
that 



J j r r 2 dfl 



dCl 



Jin 



(41) 



The optical theorem in free space relates the forward 
scattering amplitude to the cross section. 

Here, we modify the above free space relations, in order 
to define a cross section for our confined geometry. With 
incident mode <fi n (x, y), the scattering wavefunction is 



ipn{x,y) - <j) n {x,y) 




Xn(ya)xm(yo) 



X — Xm(|/) 



(42) 



which yields a probability current in the x direction 



Jx 



\xn(y')\ 2 



1 



--Xm(yo)x P (yo) 



^-^\ Xm (y) Xp (y). (43) 



Integrated over y, the cross terms cancel, and this current 
yields a scattered flux 



jxdy = I<M 

o k { x n) 



l Xm(j/0 



(44) 



going in each direction. As our basis is normalized to 
unit flux, in our case, the incident flux is simply 1. Wc 
define a differential cross section for the n th incoming 
mode, similar to the free space differential cross section 
©J: 



= dy 



dy 



\TZsVd 



da n 
dy 

J jxdy 



X. n 



„( m ) 



(45) 
(46) 
(47) 



where we note the exact analogy between plane waves in 
our waveguide and partial waves in free space. Note also 
that o n has the correct dimensions of length. 

In free space, the cross section of a spherically sym- 
metric scatterer is independent of the direction of the 
incident plane wave. In the waveguide, however, the 
cross section depends on the incoming mode, because 
the waveguide breaks spherical symmetry. We can de- 
fine a total cross section as the average over incoming 
directions, 



1 

2N 



N 



(48) 



--N 



We note that a represents a sort of fraction of the in- 
coming wavefunctions which scatter. Flux conservation 
imposes the bound 



< E a n<d. 



(49) 



n=-N 



The maximal value of a will, therefore, be 1/2N, and its 
minimal value will be 0. 

We had previously mentioned that Rank(T — I) = 
Rank(R) =1. This reduced rank indicates that, as in free 
space, a basis exists in which a single wavefunction scat- 
ters, while the remaining 2N — 1 do not. The quantity 
2No thus tells us the fraction of the flux of this scattering 
wavefunction, and should be between and 1. Hence- 
forth we will speak of a = 2Na/d as the cross section (as 
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a fraction of the wire width) . Our final cross section is 
thus 



which is very similar to the 2-D free space optical theo- 
rem, 



JV 



xlivo) 



* = w £ ,<„> 

\n=-N ft z 



and satisfies 



< a < 1. 



(50) 



(51) 



We note that for m impurities in the wire, assuming 
there are N > m modes available, Rank(T — I) = 
Rank(R) —m. By arguments analogous to the single scat- 
terer case, with m scatterers present, we could choose a 
basis where only m wavefunctions would scatter, while 
the remaining N — m would be transmitted without any 
perturbation. 



I "(«^" /4 ) = V^ 



(58) 



where /(0) is the forward scattering amplitude. Combin- 
ing (|5U|l and allows us to define a cross section, 



a(k, a) 



E 



n k x 



\Ks\ 2 



(60) 



4- Relationship between Cross Section and Conductance 



3. Modified Optical Theorem in the Waveguide 



The usual free space optical theorem involves the imag- 
inary part of the forward scattering amplitude. One can 
show that the free space optical theorem constrains s via 



(Ims) 2 = -- |s| 2 . 



(52) 



A similar relation holds in the wire, and as it turns out, 
flux conservation implies that not every value of IZs is 
physically permissible. Unitarity of S implies that 



i. 



(53) 



Using I|37I38[I in after some algebra, we find 

EJ i ?mn| 2 + <Wl + 2 Re Rmn) = 1, any n (54) 
which requires that TZs satisfy the constraint 



i^i 2 E 



Xm(Vo) 



(m) 



= -ha.(Hs). 



(55) 



We can express Ij55(l as an optical theorem, like the one 
in free space. Define the "forward scattering amplitude" 
to be the amplitude which scatters in the direction of the 
incoming wave: 



_ _iJZs_ 

In — ( n ) XnyUO) 



(56) 



With this definition, our "optical theorem" in the wire 
becomes 



(57) 



We note that our cross section is trivially related to the 
conductance of the wire, calculated as in the Landauer 
formalism^! as G = ^Tr T^T. Noting that T = I-R, 
and letting N equal the number of open channels, 



G cx TrT f T. 

= 7V + 2ReTrR + ^^|i? 



mp | 



m p 
2( 



p k x 

+ E 7M x l(yo) E T^^miyo) 

fox 



(61) 
(62) 



,2/ 



= N-cr 



(63) 
(64) 



where we have used <|3TJ) . I|55p. and JSSJ. 

In this section, we have laid out the formalism which we 
shall use to examine scattering in the wire. In Sect ion llVl 
we use our formalism to examine how the cross section 
and conductance of the confined scatterer vary as the 
physical parameters of the wire change. 



IV. RESULTS AND DISCUSSION: 
CONFINEMENT- INDUCED SCATTERING 
PHENOMENA 

Confinement-induced effects on transport, in partic- 
ular "conductance dips" as each new mode opens, have 
been observed experimentally^* as well as in other theoret- 
ical investigations^*2*iii of impurity scattering in quan- 
tum wires. Changes in the transport due to motion of a 
single defect are studied in relation to universal conduc- 
tance fluctuations^ Theoretically, the observed phenom- 
ena are generally explained via mode-mixing effects^ or 
loosely attributed to multiple scattering from the walls. 3 
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In this section, we discuss some confinement induced phe- 
nomena as we specifically observe them in our wire. 

The first quantity we wish to examine is how the cross 
section of the confined scatterer relates to the free space 
cross section of the identical scatterer, as a function of 
energy. The cross section of a free space scatterer is 

*f = {\s\ 2 - (65) 











1 














T 








1 
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Figure 4: Free space cross section (thick line) and cross section 
for confined scatterer (thin line), vs. kd. The scatterer is in 
the center of the wire, at yo = Q.5d. Its scattering length is 
a = O.ld. 

Fig. 01 shows that one feature resulting purely from 
confinement is the appearance of sharp discontinuities in 
the cross section as new modes become available; the ef- 
fect is present only if the newly opened mode is nonzero 
at the scatterer. To the left of each mode opening, the 
scatterer appears entirely transparent, whereas to the 
right of the mode opening its cross section is the full 
width of the wire. Another point to note is that, while 
the free space cross section decreases monotonically with 
k, the cross section of the confined impurity does not, 
and can be nonzero at arbitrarily high k. 

A. Resonances 

We wish to examine some of the features of the con- 
fined cross section in Fig. 0|in greater detail. In partic- 
ular, we want to look at the resonances, and understand 
the effect of varying the scatterer position. Fig. [S] shows 
how the cross section and conductance vary as functions 
of energy, for three different values of yo. We observe the 
general property that resonances in the cross section ap- 
pear as new transverse modes open (kd = mi), unless the 
newly opened mode is zero at the scatterer. We further 
observe that the resonances have a universal structure: 
If a resonance exists at kd = rm, then 

lim a(k) = (66) 

kd— m7r~ 

lim cr(fc) = 1. (67) 

kd — >?i7r+ 



As the scatterer position changes, the character, shape, 
and number of resonances change also. In Fig. the 
scatterer is located at yo = Q.05d, where each of the 
modes is nonzero. Each resonance obeys H66I67H . In 
Fig. |5J;), the scatterer position is yo = 0.32d Because 
X3(0.32d) <C 1, the resonance at kd = 37r is much nar- 
rower than the others, and will in fact vanish if we move 
the scatterer slightly to yo = d/3. In FigJSJ;), we observe 
such missing resonances-with the scatterer in the precise 
center of the wire, all even modes are zero there, and the 
resonances at kd = 2n, 4-7T, 6ir, . . . have vanished. 

Figs. Eb);Eli)i and 03) show the conductance for each 
scatterer position, which is related to the cross section by 
H64JI . Discontinuities such as H66I67J1 in the cross section 
lead to a continuously varying conductance. The impu- 
rity reduces the conductance by at most one stairstep. 
Comparing the conductance with (heavy line) and with- 
out the impurity (thin line) , we observe that to the right 
of each mode opening, the conductance falls, generally to 
the value of the next-lowest stairstep. 

From Fig. 03 our observations are as follows: (1) Reso- 
nances occur for values of k at which a new mode opens, 
unless the new mode is zero at the scatterer. (2) The 
general structure of a resonance is that the scatterer be- 
comes transparent immediately before a mode opens, and 
that its cross section jumps to one immediately after the 
mode has opened. (3) The shape and width of the res- 
onances depend on the scatterer position, yo- In Section 
fVl we shall prove these statements, and explain them in 
terms of interference effects. 

Although we have examined the case of a repulsive 
scatterer in Fig. [SJ our analytic and numerical results 
indicate that even for an attractive impurity (a < 0) the 
conductance is reduced by a single unit immediately af- 
ter the subband opening kd = n7r, rather than just below 
it. Similar results about the conductance reduction near 
the opening of each subband appear in other theoretical 
investigations of quantum wires; we cite some represen- 
tative worksASiLi& One result is for a scatterer in a 
rectangle, and therefore cannot be exactly compared to 
ours. However, the remaining referencesSi^ treat the 
infinite wire, and show a slight difference from ours: In 
Refs. 6 and 7, which examine only attractive impurities, 
the reduction in conductance appears slightly below the 
subband opening. Ref. 10 examines both attractive and 
repulsive scatterers. For repulsive scatterers, Ref. 10 
presents results similar to our Fig. [5] while for attrac- 
tive scatterers the results of Ref. 10 resemble those in 
Refs. 6 and 7. Because delta functions in more than one 
dimension do not scatter, many different methods exist 
of representing a point impurity in two dimensions. We 
have chosen to use a t matrix representation as in Refs. 
12, 16, and 17, and the particular t matrix we have cho- 
sen in mimics the s wave scattering of a hard disk of 
radius a. The references whose results differ from ours 
have used explicit limiting forms of delta potentials. We 
attribute the slight difference in results to the different 
forms and representations of the point impurity. 
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Figure 5: Cross sections (a,c,e) and conductances (b,d,f) vs. the wavenumber, k, for the scatterer at three different positions 
in the wire. The conductance is in units of the conductance quantum, 2e 2 /H. The scatterer positions are yo = 0.05d (a, b), 
j/o = 0.32d (c, d), and yo — 0.50d (e, f). The scattering length is a — O.ld. This figure illustrates resonances, and how they 
change with scatterer position. Resonances in the cross section appear as new transverse modes open (kd — nn). In general, 
just below each mode opening, a — 0. Just above each mode opening, <x=l. By moving the scatterer around, we can tune the 
character of the resonances. In (a), the scatterer is at yo — 0.05. Since Xn(0-05) 7^ for n any available mode, each resonance 
obeys the general rules 1661671 . In (b), because X3(0-32d) <C 1, the kd — 3tt resonance is very narrow, and barely there. In (c), 
with the scatterer in the precise center of the wire, Xn(0-05d) 7^ where n is any even available mode, and thus the resonances 
at kd = 0, 2-7T, 4-7T, . . . have vanished entirely. Discontinuities in the cross section lead to a continuously varying conductance: the 
two quantities are related by (I64H . The impurity reduces the conductance by at most one channel. Comparing the conductance 
with (heavy line) and without the impurity (thin line), we observe that to the right of each mode opening, the conductance falls, 
generally to the value of the next-lowest stairstep. Several referencea2iS*£*— obtain slightly different results for the conductance 
vs. wavenumber curve, in which a — 1 just below the mode opening. We attribute the discrepancy to differences in how the 
point scatterer is modeled. 
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B. Between Resonances 

We have observed that the scatterer position in the 
wire influences not only the width and shape of the res- 
onances, but can even cause a resonance to disappear. 
Between resonances, we note the opposite phenomenon: 
The scatterer position barely affects the cross section at 
all. As we vary the scattering length a, the cross section 
behaves very similarly to the free space cross section (see 
Fig. 0. 



that in the waveguide, as in free space, a basis exists 
in which only a single wavefunction scatters. We write 
this basis explicitly, and show that the scattering of this 
single wavefunction allows us to explain the confinement- 
induced phenomena shown in Figs. 1 ll(il) simply, in terms 
of reflections from the scatterer and its images. 



A. Single Scattering Wavefunction and the Hall of 
Mirrors 




-0.1 -0.05 



Figure 6: (a) Cross section of the confined scatterer for 
kd = 12.57T (halfway between having 12 and 13 modes open). 
On the horizontal axis, we vary the effective scattering radius 
a from — 0. lei to O.ld. The vertical axis is yo, the vertical posi- 
tion of the scatterer in the waveguide: yo varies from the edge 
of the wire (yo = 0) to the middle of the wire (yo = 0.5). In 
contrast to the case of a resonance, between resonances, the 
cross section is almost identical to the free space cross section 
af, shown in (b); the arrow indicates two corresponding max- 
ima. The only effect of varying yo is to introduce some small 
oscillations due to interference effects. The plots continue in 
a similar manner to the left and right. 

Our general observations from the numerical results 
are the following: (1) Scattering resonances occur where 
nonzero modes open. (2) The cross section jumps from 
zero to d across a resonance. (3) The width and exis- 
tence of a resonance is influenced by tjq. (4) Away from 
resonances, the cross section behaves like the free space 
cross section. 



V. PHYSICAL PHENOMENA EXPLAINED IN 
TERMS OF A SINGLE SCATTERING WAVE 
FUNCTION 

We have shown that items 1 1121 in our list of Section 
III Bl are characteristics of free space scattering which have 
analogs in the wire. We have yet to examine whether 
item |3 has an analog as well. In this section, we show 



In the basis of plane waves, the scatterer couples the 
channels so that, in general, an incoming mode scatters 
into all the modes. However, 



rank 



S - 



I 

1 



= 1, 



(68) 



which implies that a particular choice of basis exists in 
which only one of the incoming wavefunctions scatters at 
all, and is decoupled from the other N— 1 basis functions. 
This wavefunction is the analog of the s wave in free 
space. 

Another way of seeing that only one wavefunction scat- 
ters is by the following argument: For a given energy, N 
channels are open. Any incoming wavefunction is a linear 
combination of the N basis functions 



4> {n \x,y) = 



-Xn(y), 



(69) 



(we require only the left- or right-moving set). Any linear 
combination 



N 



(70) 



71=1 



of the N basis functions Ij69(l which is nonzero at the 
scatterer will satisfy 



N 



X)cn0 Cn) (3O,!to)^O. 



(71) 



Equivalcntly, define an 1 x N matrix W such that W n = 
4>( n \%o,yo) and a vector of coefficients c = [c\, . . . cjv]- 
Then, for any scattering wavefunction, the c n will satisfy 



WcV 



(72) 



As W is a rank 1 matrix, only one solution exists. That 
solution is simply c£ Nul(W) 1 - = Row(W) = W*. That 
is, the single scattering wavefunction is 



N 

<j>f(x,y) = £(^ (n) (*o,2/o)) 4 in \x,y) 



(73) 
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N 1 

= ~(n)Xn(yo)Xn(y)e ±lk 



n=l ft z 



< n) (*-*o) 



(74) 



where the choice of sign determines whether the incoming 
wave is right or left-moving. Suppose we take a linear 
combination of the two so the scattering wavefunction is 
symmetric in x. The scattering wavefunction becomes 



N 



1 



<t> s {x,y) = ^2— Xn{y)xn(yo) cos k x n \x-x ) (75) 

n=i k x 

(the other linear combination doesn't scatter). By com- 
parison with we recognize the above form of the 
scattering wavefunction as 



4> s {x,y) = -lmG w (x,y;x ,y ;k). 



(76) 



Note that all the definitions in this section, and in partic- 
ular 1)76(1 . are independent of the specific form of X n (y), 
and thus apply to arbitrary guides, not only to the hard 
wire. As we shall discuss later, $7ty is particularly im- 
portant because it shows that for an arbitrary guide, we 
can obtain the scattering wavefunction entirely from the 
Green's function, which we can approximate semiclassi- 
cally. 

We now return to the specific case of the hard wire, and 
try to understand the scattering wavefunction physically. 
From (|18J) . an alternate form of the Green's function is 

oo 

G w (f,r ,k) = - l - (-l) n H (k\f-f n \) (77) 

n=— oo 

where the f n are the image positions. Combining l|76|) 
and (|77|l . we find an alternate expression for the scatter- 
ing wavefunction: 



1 



OO 



4>s{x,y) = - ^ 

71— — OO 



(-l) n Jo(k\r-f \). (78) 



Eq. H78fl shows that the analog of the free space s wave 
in the hard wire is, as we might have expected, simply 
the free space s wave, plus an infinite series of images-a 
"hall of mirrors s wave." We can complete the analogy 
between free space and our wire by expressing the N — 1 
unscattered wavef unctions in terms of a "mirrors" basis. 
We do so by including higher order mirrored waves as in 
Tabled 

The mirroring, and the signs assigned to each image, 
are somewhat subtle: Only those combinations of free 
space orbitals which satisfy the boundary conditions, but 
are also zero at the scatterer, will not scatter. We note 
that the basis in Tabled is not an orthogonal one; inner 
products of an p wave centered on one image and a d wave 
centered on another image, for example, do not vanish, 
nor do the various such cross-terms cancel each other. 
However the important point is that the higher partial 
waves we have defined are all orthogonal to the mirror 



s wave: Higher mirrored partial waves do not scatter, 
whereas the mirror s wave does. 

We note that the image sums in Table[I]converge slowly 
and are unsuitable for numerical purposes. However, we 
have the plane wave expansion (|75|l for the s wave, as 
well as the analogous expansion 



1 



tiW = \ £ UW~r n \) (70) 

n— — oo 
n N i 

V 1 fnny\ fmry \ 



=1 k 



X COS 



fcjW(ss-so) 



(80) 



for the s wave plus positive images. We can apply the 
raising operator 



La 



1 



(d x + id y ) 



(81) 



to these plane wave expansions to obtain plane wave ex- 
pansions for the higher partial waves. For example, we 
can obtain plane wave expansions of the waves in Table 
HI as 



P Ar) = Re 



L+(i> s {r) 
L 2 +<j>f(r) 



\d x (j) s (r) 



2 >d xv (i>t(r 



k 2 ^xyYs 



and we could proceed similarly to obtain higher partial 
waves. 



B. Interference Effects 

We might expect that the cross section is related to 
the fraction of the single scattering wavefunction which 
actually scatters, and we might thus expect the cross 
section to be proportional to the amplitude of the mirror 
wavefunction at the scatterer. However, some numerics 
reveal that the mirror wavefunction alone is insufficient 
to fully describe the scattering. From (|5U|I and l|75|) . we 
find 



|^s| 2 ^(f ) 2 



s<t>s(r ) 



(82) 
(83) 



meaning that the relevant quantity is the renormalized 
mirror wavefunction, and that its value at the scatterer 
alone explains all the scattering phenomena. Eq. I|83|l is 
again independent of the form of the waveguide. 

We have already discussed how renormalization of the 
scattering wavefunction arises from interference between 
all possible bounces off images. The implication of (|83|l is 
that the value of the resulting interference pattern at the 
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Mirrored wave 


Expression 


Wavefunction plotted for kd = 40 (twelve modes open) 
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Table I: Mirrors basis. The first column indicates the symmetry, the second column the expression, and the third column a 
plot, of each mirrored wave. The arrows in the plots indicate the impurity, located at yo = 0.6d. The mirror s wave (first row) 
is nonzero at the impurity position f*o, and thus scatters. The scatterer is, however, transparent to the higher partial waves 
shown in the second through fourth rows. These partial waves, expressed as p x ,d xy ,f x 3_ 3xy 2 waves plus images, do not scatter 
due to their nodal lines along x — xq (indicated by the dashed white lines in the plots). 



location of the scatterer determines the cross section: At 
configurations (wire widths, incident wavenumbers, etc.) 
for which this interference is destructive at the scatterer, 
the scatterer is transparent. If the system is configured so 
that the bounces interfere constructively at the scatterer, 
a resonance results. 



Semiclassical Approximations for General Wires 



We wish to highlight another important point, which 
is that combining (|83|) with and (JBJ, we find that 




sJmG w (r,r ) 



1 - s [G w (r, f ) - G (r, f )} 



(84) 



The significance of (|84p. which applies to an arbitrary 
wire, is that the cross section is fully determined by 
the Green's function of the empty guide. This relation 
is important because the Green's function is a quantity 
which we can approximate semiclassically. Substituting 
the semiclassical Green's function into (|84|l will give us a 
semiclassical approximation to the cross section. Semi- 
classical methods may be of use in determining the trans- 
port properties of guiding potentials which are more com- 
plicated than the hard wire, and thus not amenable to 
exact quantum treatment. 



Scattering Resonances and Conductance 
Reduction 



In Figs. 1 4151 we observed a reduction in the con- 
ductance, which are similar (although not identical) 
to the conductance dips observed in other theoretical 
investigations^SiLiS Using our formalism, we can under- 
stand the resonances and conductance dips of Figs. 14151 
simply, in terms of interference between different paths. 
Again the results of this section do not depend on the 
specific form of the transverse modes Xn(y), and conse- 
quently apply to a general waveguide. 

As we showed in (|32|l . the renormalized mirror wave- 
function contains contributions that have propagated 
from each of the image scatterers. We are interested in 
the behavior of a cross section when a mode is about to 
open. The mirror wavefunction involves only open chan- 
nels, and is thus discontinuous across a mode opening. 
Using (|75|l . we can examine the behavior of the mirror 
wavefunction near kd — Ntt, where the N th mode opens. 
Suppose that the newly opened mode is nonzero at the 
scatterer, xn(vo) ^ 0. On the left of the mode opening, 



lim^s r ;k 



Ntt - e 



d N -' ~* 



y (85) 
7T ^ y/N 2 - n 2 



is finite. Immediately after the mode opening, 



Ntt + e 

hm (j) s [r ;k = 



d x 2 N (yo) 

7T V2A7 



(86) 
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diverges as e" 1 / 2 . That is, the mirror wavefunction is 
always finite immediately before a nonzero mode opens, 
and diverges afterwards. The value of X%(yo) determines 
the width of the resonance (and the lifetime of the cor- 
responding quasibound state). 

However, while the limiting behavior of the mirror 
wavefunction influences the shape of the resonances, it 
does not fully describe their shape. As we discussed, the 
general structure of a resonance is that a drops to zero 
just before the mode opens; see e.g. Fig. CD Clearly 
the incoming wavefunction, which includes only modes 
which are already open, cannot explain this transparency. 
The renormalization factor due to the wire, however, in- 
cludes all modes, both evanescent and propagating-and 
unlike the incoming mirror wavefunction, varies contin- 
uously across the mode opening. Using the expression 
(|B25|) for G r , we find that 



lim G r 



lim G r 



k = 



Nn-e 



k = 



Nn 



d x 2 N {yo) 

id x 2 N (yo) 

7T ^/2N~E 



(87) 
(88) 



generally. Combining 1|28|) with the limits I|85I88[) . and 
using JH2J), we can show that 



lim a \ k = e | = 

+o V d 



2Ne 



y-^ 1 [Xn(ya)/XN(yo)} 2 



(89) 



The value of Xn{uo) controls the width and structure of 
the resonance. Examining the right side of the resonance, 



lim a I k = — — 

e^o \ d 





s 


) = lim 


1 - sG r 





4> 2 s (r ) 

\Gr\ 2 

1 



^(r )(90) 



(91) 



under the assumption that \sG r \ S> 1. An important 
point to note is that when \G r \ — > oo, as it does at a 
resonance, the free space properties of the scatterer have 
little effect on the cross section (see for example 
and so the structure of resonances is universal. 



Phase Shift and Ramsauer-Townsend Effect 



We would like to note, in passing, that the trans- 
parency of the scatterer at certain energies is reminiscent 
of the Ramsauer-Townsend effect in free space. In free 
space, we define the s wave phase shift do of an s wave 



scatterer by 



1>(f) = 4>(ro) + e% \ - g (1) (fc \r~ n>|), (92) 



Let us define it in the confined scatterer as the phase 
shift of each s wave generated from the source and the 
images, when the incident wavefunction is our mirror s 
wave: 



e 2i(5 _ ^ 



(93) 



x (fc|r--f n |) 



Comparing to the form (|26|l of the Lippmann-Schwinger 
equation, we find that 



,2i<5 



1 - is<f) s (ro). 



Applying 1)83(1 . we find that 



,2i<5 



(94) 



(95) 



so that a = when <5o = 0, n and a — I when So = § , 
As in the free space Ramsauer Townsend effect, the cross 
section vanishes when our analog of the s wave phase 
shift does. 

In conclusion, we can make the following general state- 
ments about the cross section: 

1. If Xn{vo) 7^ 0, the cross section will be discontin- 
uous at kd — Ntt. The limit on the left hand side 
will be finite, and depend the value of xn(vo)- The 
limit on the right hand side will be unity. 

2. The width and shape of the resonances in item (2) 
will depend on Xn(,Uo)- Physically this means that 
one can tune the resonances by sliding the scatterer 
up and down in the wire, or by changing the nature 
of the confining potential. 



D. Semiclassical Interpretation of Resonances 

Having explored the structure of the resonances, we 
wish to understand them semiclassically, from an inter- 
ference point of view. We again turn to the hard wire for 
insight into the scattering processes. The scattering cross 
section depends only on a single quantity, the renormal- 
ized mirror wavefunction. We have previously expressed 
this renormalized mirror wavefunction as a sum over the 
different paths ending on the scatterer l|32() . When the 
interference between these paths is maximally construc- 
tive, the cross section is maximal. When it is maximally 
destructive, the cross section vanishes. 

In order to make this interference clearer, we shall 
make a simple approximation. Consider the most 
straightforward case, with the scatterer in the center of 
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the wire (the case we examined in Fig. |2J. The image 
positions are r n = (0,nd). In the Green's function lfT5)l. 
we replace each Hankel function by its asymptotic form, 
which is equivalent to making a semiclassical approxima- 
tion. This yields an approximate Green's function 



G w (r,f ) = 



1 



= 5ttz/4 



27T 



E 



1 



n— — oo 
i(k\r—r n | — 7i.7r) 



(96) 



The contribution from each image carries a phase related 
to the path length from the image to the observation 
point, as well as a Maslov index e m7r describing the sign 
change after n reflections from the wire walls. In the limit 
that r — > r*o, we see that the contribution propagating 
to the scatterer from the n th image has relative phase 



i{kd 7r)n p or yames f sucn faat 



kd = (2p — 1)tt, p integer, 



(97) 



all the scattered wavelets thus interfere constructively. 
Comparing with Fig. 0] we see that this constructive 
interference coincides precisely with the resonances which 
occur as new modes open. 

To fully explain the structure of the resonances, 
and understand the scatterer's transparency just before 
modes open, we would have to consider the effects of 
renormalization in this approximation. We shall not pur- 
sue the semiclassical limit further here, as we have al- 
ready examined the exact case in great detail in Section 
IV CI However, purely from this simple semiclassical ar- 
gument, we can see that the resonances induced by the 
wire are indeed related to interference effects between 
different reflections. 



VI. CONCLUSIONS AND FUTURE 
DIRECTIONS 

Combining the method of images with the idea of a sin- 
gle scattering wavefunction, we have developed a formal- 
ism with which to treat scattering from a single impurity 
confined in a quantum guide. We find many similarities 
between scattering in the confined geometry and scatter- 
ing in free space: We have defined meaningful analogs 
of the free space cross section, optical theorem, partial 
waves, and Ramsauer-Townsend effect. Although the 
hard guide is particularly useful for insight, many of our 
results apply to other types of waveguides also. 

We have examined the transport properties of a 
confined scatterer, making several general observations 
about the existence, locations, and and characters of the 
resonances. Additionally, we have outlined a method for 
approximating the cross section and conductance semi- 



classically, for arbitrary guides. We have derived an effec- 
tive optical theorem for general quantum wire potentials. 
Using the fact that a confined target is rank one, we have 
described the transport properties of the wire entirely in 
terms of a single scattering wavefunction, renormalized 
by reflections from the confining potential. Our central 
result is that the cross section of a confined impurity is 



S(f> s (r ) 



where 4> s (r) is the renormalized single scattering wave- 
function, which includes interference effects due to reflec- 
tions from the waveguide walls. 

Despite the common use of the hard-walled guide as 
a theoretical model for 2DEG quantum wires 
actual confining potentials for both atom waveguides as 
well as 2DEG quantum wires tend to be soft or even 
parabolic-whereas in nanotube quantum wires, the rel- 
evant boundary conditions are periodic (the periodic 
boundary case is also exactly solvable via the method 
of images, although we have not presented the deriva- 
tion here). While the hard wall itself is not the best 
model for 2DEGs, it is the ideal system in which to ex- 
amine how effects like conductance reduction, which are 
independent of the confining potential and occur in more 
realistic geometries, arise from simple reflection and in- 
terference phenomena. The effects of multiple scattering 
from the impurity are in fact very similar for general 
confining potentials (see e.g. the renormalized t matrix 
formalism^) . 

A second limitation of this work includes the inade- 
quacy of the s wave scatterer model at high energies, and 
in particular at energies which are sometimes relevant 
when imaging electron flow. In a future publication, 22 
we plan to extend our formalism to treat higher partial 
waves. Another possible extension of this work could be 
to examine the transport properties of many scatterers in 
the wire; for example, explicitly explaining the transition 
from ballistic to diffusive transport in terms of interfer- 
ence between multiple mirror wavefunctions. 
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Due to the linearity of Schrodinger's equation, we can 
consider each of the constituent plane waves as scatter- 
ing independently. The relevant physics thus reduces to 
diffraction of a single plane wave incident on a periodic 
array of scatterers. 

When our incident wavefunction strikes the diffraction 
grating, we expect that in the Fraunhofcr limit, far from 
the array, the diffracted beam will be a sum of plane 
waves at the Bragg angles (see Fig. 0). 



A 




Figure 7: A mode incident on an impurity in the hard wire is, 
effectively, a superposition of two plane waves incident on an 
infinite diffraction grating. At large distances from the scat- 
terer, the scattered wave must therefore be a sum of plane 
waves at the Bragg angles. Closer to the scatterers, evanes- 
cent modes appear also. 



Appendix A: WIRE AS A DIFFRACTION 
GRATING 

In Sec. IIII Al we used the mathematical equivalence of 
the wire to a periodic array of image scatterers. Further 
extending the analogy, we now reformulate the problem 
of the confined impurity as one of diffraction from an in- 
finite, periodic grating, and the scattered wavefunction 
as a diffracted beam. We wish to note the detail that be- 
cause of our Dirichlet boundary conditions, our incident 
beam is a sum of two plane waves rather than a single 
plane wave: 

4> m {x,y) = —^e ik ^- x °\ X m{y) (Al) 



Each order of the diffracted beam will have a different 
weighting; our object is to cast our Green's function ltTH|) 
in a manner that makes these weightings explicit. This 
formulation is of interest because it is equivalent to de- 
termining the transmission and reflection coefficients of 
our system. Namely, if we send in a mode, we would like 
to know what modes come out and with what intensi- 
ties, which is precisely what we shall find by casting our 
scattering wavefunction as a diffracted beam. 

To begin, we note that Bessel functions are regular, 
and consequently have plane wave expansions. In fact, 
the well-known integral form of the Bessel function 

Jo(kr) = ^ '^ m ' r 'de (A3) 
27T Jo 

expresses the Bessel function as an isotropic sum of plane 
waves emanating in all directions from the origin. We can 
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rewrite (|A3f) in terms of a line integral over k y a.s 



Jo(kr) 



1 

1 

2^ 



d(ksmd) Hk x ,k v ).{x,v) 



kcosO 



dk 



(A4) 
(A5) 



where C lies along the real axis, going back and forth 
within the interval (— k, k). 

We seek an analogous expression for the Hankel func- 
tion. Since \k y \ is always less than k in i|A5[l . k x is always 
real, and the plane waves in I|A5J| do not include evanes- 
cent waves. Suppose we modify the right-hand side of 
(|A5I) . By permitting k y to range from — oo to oo, we 
may include contributions from all possible evanescent 
waves, yielding 



1 

2^ 



dk 



(A6) 



For 

where k x 



k y \ > k, the corresponding value of k x — ±in x , 



1.2 

K v 



k 2 , is imaginary. Choosing the posi- 
tive value makes the evanescent waves die off to the right 
(x > 0) and diverge to the left [x < 0). Therefore, with 
this choice of sign, we may expect (|A6ll to converge only 
in the right half-plane. 

By adding on an evanescent contribution to the Bessel 
function's plane wave expansion, we have in fact arrived 
(up to a constant) at an integral form of the Hankel func- 
tion, valid in the right half-plane only, 



where 



H^\kr) 



k x — 



dk 

k x 



h' 2 

K v 



\ky\ < k 



l\/hy hP 1 |^y| 



(A7) 



(A8) 



Eq. (|A7|> is a well-known expansion, 20 readily obtained 
by Fourier transforming Green's equation. 
Consider now the sum 



G P {r-d) = - 1 - J2 H (k\r-(x ,nd)\) 



(A9) 



n— — oo 



which is the Green's function for a fully periodic array 
of scatterers at (xo,nd) (or equivalently, a periodic wire 
with walls at y = 0, y = 2d and a scatterer halfway be- 
tween). Based on the physical picture of diffraction (Fig. 
0), the sum in (|A9|I must be a superposition of plane and 
evanescent waves at real and complex Bragg angles re- 
spectively. We wish to find the weighting on each wave 
(analogous to the structure factor in X-ray diffraction). 

To find the weights, we calculate the expansion explic- 
itly. Using in (|A9|) . we find 



G p (r,0) 



E H ^ 



n— — oo 



I 

2^ 



E 



n— — oo 



dk y 

~k~ 

— OO "'X 



xe 



i(k x ,k y )-(x — xo,y—nd) 



I 

'2?r 



dk 



k 

— OO "'X 



y_ e i(k x ,k y )-(x-x ,y) 



(A10) 
(All) 

(A12) 



OO 

X 

n— — oo 

dk 



oo 

j2 (e-wv^y 



k 

-oo ">x 



y_pi(k : c,ky)-(x-x ,y) 



(A13) 



2m7i\ 



m— — oo 



where we have used 

oo oo 

E [e lfc « d ]™ = 2tt E S(kyd-2im), 



(A14) 



n— — oo 



n— — oo 



(similar to stationary phase) to reach the expression 
(|A13(1 . One can, more rigorously, reach l|A13() via the 
Poisson sum formulai22i2i 

Performing the integral in 1A13J1 yields the final ex- 
pression 



G p (r, r m ; k, d) 



oo 

* y J_ e «i n) l*-»ol(A15) 
d ^ u( n ) 

-—oc ^x 



n— — oo 



x cos(fc( n) y) 



where we use the absolute value signs to extend the 
Green's function to converge on x < xq, and we define 



k x n) = fccos0„ 
k y n) = ksin9 n 

where the 9 n are the Bragg angles 

2mr 



= arcsm ■ 



kd 



n integer 



(A16) 
(A17) 



(A18) 



and the complex- valued 6 n are the result of the logarith- 
mic singularity. 

The image representation of our wire is not exactly the 
periodic array of Fig. [7] The images alternate in sign, 
and in general the scatterer is off-center. We can, how- 
ever, represent our image array as a sum of two arrays, as 
in Fig. [S| One can easily verify that the Green's function 
for the wire with point source at = (xo,yo) becomes 
the sum 



G w (r,r ) 



[G p (f- y y;2d) 



(A19) 
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2d 


'o~+ 





O - 


2d 






'o - 



where a is a constant which we shall determine. The 
static Green's function happens to be the potential which 
solves Poisson's equation for a point charge between two 
grounded, conducting plates^ and we can put it in 
closed form as follows: 

OO j 



Figure 8: In the case of the off-center scatterer, the relevant 
image array is the sum of two periodic arrays of opposite 
signs. If the source is at (xo,yo), the image locations are 
(xo,yo + 2nd), (xo, —yo + 2nd) for the positive and negative 
arrays respectively, where n — ±1, ±2, . . .. 



-G p (f+y y;2d)] 

OO 

^ An) e 

n=l h x 

x Xn(y)xn(yo) 



(A20) 



where the Bragg angles are modified slightly from the 
case of the periodic wire, and defined by 



smf,, — 



kd 



(A21) 



The form l|A20(l is identical to (|14fl . We have shown that 
diffraction translates the method of images into the usual 
eigenfunction expansion of the Green's function. 



where 



1 d 2 . /mny\ 
= — > — x - sin — — 
7T * — ' ma V a / 



(B4) 



1 °° 1 

1 E 1 

rn—l 



d J 
'mw(y + y ) 



cos 



1 °° 1 
= —Re — e - ^ 



(B5) 

e -^|a;-a;o 

(B6) 



m— 1 
X |f" 
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T (y + yo) im7T(y-y ) 

d — e a 



Z± = e ^ x x °'e 3 



Using the identities 



1 / \ 

= iRe E ( — - — ) ( B7 ) 

7T A — ' V 771 ml 

TCI — 1 V ' 



(B8) 



Appendix B: IMPROVING CONVERGENCE OF 
THE GREEN'S FUNCTION 



E 

in. — i 



ln(l - Z) 



(B9) 



Form l|14l) of the Green's function converges slowly 
also. The reason for the slow convergence is the sin- 
gularity at the source, which is present but disguised in 
the evanescent modes. We use the Kummer method for 
convergence acceleration^ which also makes the singu- 
larity explicit. Kummer's method involves adding and 
subtracting a multiple of the static (k = 0) Green's func- 
tion. We note that [12] obtains similar results, via a 
slightly different procedure. 

The wire Green's function is 



G w (f,r;k) = -i ^ -^Xm(y)Xm(yo) (Bl) 

m=l 

o ik^"-'\x-x \ 



xe 



We wish to consider a more rapidly converging form of 
the Green's function: 

G w (f,f ;k) = [G w (r,r Q - 1 k)~aG w (r,r o - 1 0)] (B2) 
+aG w (f, f ;0) 



and 



we find 



RelnZ = hx\Z\ 



G w (f,f ;0) = -In 

7T 



1 - Z_ 



In 



2vr 



= 7^ In 
2tt 



= — In 

2tt 
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- z_ 
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-z+ 





\x— xo\ 



e f \x-x \ _ e ^-(y+yo) 
e f \x-x \ _ e f (y-yo) ' 



\x—x \ 



e d 



^(y+yo) 



\x-x \ _ g-^ty-yo; 



\x—xo\ 



e a 



(y+yo) 



(BIO) 

(Bll) 
(B12) 

(B13) 
(B14) 
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1 , cos 
in 



[f(l/- yo)] - cosh [f(a:-xo)] 



2tt cos [§ (y + ?/o)] - cosh (a: - x )] ' 

The singularity of the Green's function arises from a 
logarithmic singularity in the Hankel term at the source. 
Near the scatterer, the contribution from this term be- 
comes 



lim Go(r,r ; k) = 



r \ (B16) 



lim H (k\f-fo\) (B15) 

1 

+ -7 

7T 

7 - In 2 



(B17) 



In (k\r- r |) 



where 7 is the Euler-Mascheroni constant, while from Eq. 
(|B14|) . wc find that the limiting behavior of the static 
Green's function is (after some algebra) 



Jim G w (f,f ;0) 



--In — sin C^f) (B18) 
-\ — In (k\r - r \) . 



In order to cancel the logarithmic singularity in Green's 
function, we thus require that 



a = 1 



(B19) 



Our final form for the Green's function, using (|B2IB19| > 
is 



G w (f,r ; k) 
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Xm (y)x m (yo) 



1 



a ikl m) \x-x \ 
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, 1 ln cos W$ z Vo ^ z cosh [%( x ~ 

2vr cos [i(y +y )] -cosh [^(x - x )] 

This expression is of further use in that it allows us to ob- 
tain another form of the renormalization constant 1231) . 
which we had shown to have the slowly convergent ex- 
pansion 



(B22) 



G r = --J2(-V m H (k\r m -r \). 

We can now express G r in an equivalent, but more 
rapidly convergent expression, suitable for numerical pur- 
poses: 

G r = lim {G w {r, r ;k)- G (r, r ; k)) (B23) 
= lim [(G w (r,r ; k) - G w (f y f ;0)) 

r— H- 

+ (G w (r,r o ;0)-G o (r,r ;fc))] (B24) 
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so that 



lim [Gu,(f, rb; 0) - G (r, r ; fc)] 
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